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Abstract. We study the hydrodynamic interactions between colloids suspended in a 
compressible fluid inside a rigid channel. Using lattice-Boltzmann simulations and a simplified 
hydrodynamic theory, we find that the diffusive dynamics of density perturbations (sound) in the 
confined fluid give rise to particle correlations of exceptionally long spatial range and algebraic 
temporal decay. We examine the effect of these sound-mediated correlations on two-particle 
dynamics and on the collective dynamics of a quasi-one-dimensional suspension. 



1. Introduction 

Particles embedded in an unbounded fluid are dynamically correlated by hydrodynamic 
interactions that decay like 1/r (where r denotes the separation between particles) pQ. This 
result is understood in terms of the fundamental solution to the steady-state Stokes equation, 
r]X7 2 u = Vp — b(5(r), supplemented by the fluid-incompressibility condition V • u = 0. Here u is 
the flow velocity field, p the pressure field, b a point force, and rj the shear viscosity. It follows 
from the first equation that the vorticity, V x u, satisfies the Laplace equation, leading in the case 
of an unbounded fluid to the 1/r slow decay of the velocity field. When particles are confined to 
a quasi-one-dimensional (qlD) channel, however, these long-range dynamic correlations are cut 
off [2], the screening being attributed to dissipation at the fluid-solid interface, i.e., to the loss 
of fluid momentum at the channel boundaries. The resulting screening length is comparable to 
the channel width. 

The dissipative effect of the boundaries can be incorporated through an additional friction 
term in an "effective fluid" model ]3l 0], ryV 2 u = po^u + Vp — b<5(r). In the extra term po 
is the fluid mass density and £ an effective friction coefficient. The parameter ^ has units 
of inverse time and characterizes the rate of momentum loss to the boundaries. It can be 
estimated, therefore, as the inverse of the time it takes fluid momentum to diffuse to the channel 
boundary, £ ~ v/h 2 , v = r]/po being the kinematic shear viscosity and h the channel width. This 
simplified, phenomenological approach has been shown by simulations [5j El [8] and analytical 
calculations [9] to correctly reproduce the qualitative behavior of confined fluids. The modified 
Stokes equation leads to a Helmholtz equation for the vorticity, whose regular solutions are 
exponentially screened functions. Yet, the friction that suppresses the vorticity (transverse fluid 
stress) does not similarly suppress the pressure field (longitudinal fluid stress) emanating from 
a locally applied force. Consequently, the application of a localized force generates a long-range 
pressure distribution, which may give rise to long-range flows and forces on embedded particles. 
For example, in a q2D suspension confined between two plates, the pressure distribution creates 



hydrodynamic interactions that decay like 1/r 2 , while the 1/r correlations due to transverse 
momentum transfer are cut off [10]. Thus, although the hydrodynamic interactions in the q2D 
system remain long-ranged, the very mechanism for correlations changes [llj. In the case of a 
qlD channel, the pressure generated by a localized force far along the channel is constant and, 
therefore, does not produce long-range forces. Thus, within the Stokes description of a steady 
incompressible flow, in a qlD rigid channel long-range hydrodynamic interactions arise neither 
from transverse stresses nor from longitudinal ones, and the correlation between particles is 
screened, the screening length being proportional to the channel width [2J. 

The omission of compressibility, and with it of compression (longitudinal sound) modes, is 
usually justified for unbounded fluids, where the effects of sound are short-lived. The short 
propagation time to any reasonable distance r, t s ~ r/c s , is a result of the large speed of 
sound, c s ~ 10 3 m/s for ordinary liquids such as water. The effect of sound on the short-time 
hydrodynamic interactions in unconfined suspensions was addressed in earlier works [12 1 113 1 PH], 
The corresponding effects are subtle and short-lived, quickly taken over by the fully established 
hydrodynamic interaction due to transverse flow. 

In confined geometries, however, the behavior of sound modes qualitatively changes. Inclusion 
of fluid compressibility has been found to lead to a velocity autocorrelation function of a confined 
particle, which decays with time only algebraically, with a negative long-time tail [5| |6| \7\ [To]. 
[In the case of a qlD channel the decay is ~ (— i~ 3 / 2 ).] Thus, under confinement, sound 
acquires long-time memory without a well-defined relaxation time, and the incompressibility 
assumption loses, at least in principle, its justification. Informed by this result, one suspects 
that the long memory of sound modes may play a role in the dynamic correlations between 
different particles in a channel: since sound may propagate to long distances, the resulting 
correlations will be long-ranged, and the algebraic temporal decay will endow these interactions 
with a long-time memory. This novel mechanism for correlations challenges the conventional 
wisdom described above, according to which fluid-mediated correlations in a rigid channel are 
exponentially screened. The question that comes to mind is whether the long-time memory 
of sound is sufficiently significant to be observed experimentally and to undermine predictions 
of the theories based on the assumption of incompressibility. We remind that the algebraic 
temporal decay of sound modes in a channel does not imply that sound is capable of causing 
steady rearrangement of the fluid. The integrated effect of the i _3//2 tail decays as t _1//2 , which 
is very slow, but still implies that at infinite time (the steady-state limit) all perturbations due 
to sound eventually vanish. 

In a recent Letter [16] we have demonstrated that the heuristic arguments given above are 
correct — namely, the diffusive, slowed-down sound modes of a fluid confined in a rigid channel 
mediate long-range and long-time velocity correlations between suspended particles. In the 
present publication we report detailed findings concerning the compressive response of the fluid 
and its effect on particles, which were not included in that brief publication. Two tools are 
employed and compared: simplified analytical calculations and lattice-Boltzmann simulations. 
We assume no-slip boundary conditions at the fluid-particle interface and consider either slip 
or no-slip boundary conditions at the channel boundaries. Fluid dynamics in a qlD channel 
strongly depends on the boundary conditions imposed at the edges of the channel [T7]. In the 
calculations we assume an infinitely long channel with a vanishing flow velocity at the edges. In 
the simulation we take a sufficiently long channel, such that all the temporal results presented 
below are insensitive to the system size. Further details of the simulations can be found in 

Ref. ma. 



2. Dynamics of fluid density perturbations 

2.1. Unbounded fluid 

We begin by revisiting the compressive response of an unbounded Newtonian fluid to an 
impulsive force [18], to which we will later compare the case of a fluid confined in a channel 

[EJUS]. 

Isothermal fluid dynamics [20] is governed by the momentum-conservation Navier-Stokes 
equation, 

pd 4 u+(u-V)u =-Vp + 7 ? V 2 u+(7 ?1 , + r ? /3)V(V-u)+f, (1) 
and the mass-conservation equation, 

d tP + V • O) = g. (2) 

In these equations p is the fluid density field, f an external force density, and r\ and rj v are 
the shear and volume viscosities, respectively. The function g could represent, for example, a 
mass monopole, m5(r — ro), where at the point ro fluid is generated or lost, or a mass dipole, 
— d • V<5(r — ro), where fluid is at the same time created and lost, thus globally conserving mass 
but generating a flow. 

The hydrodynamic equations are linearized by neglecting the convective term, (u • V)u, 
and separating the thermodynamic variables into equilibrium and small perturbation parts, 
p(r,t) = po + Sp(r,t) and p(r,t) ~ po + (dp/dp)5p = po + c 2 5p(r,t). The resulting linearized 
equations read 

p d t u = -c 2 V5p + rjV 2 u + (rj v + r?/3)V(V • u) + f , (3) 

d t Sp + po V • u = g. (4) 

When dealing with sound, it is convenient to decompose the linearized equations into 
transverse and longitudinal components (the so called Helmholtz decomposition), u = u^ + u'", 
which satisfy V • \i T = and V x u L = 0. The decomposition yields 

p d t u T = r?V 2 u T + f T , (5) 

and 

pod t u L = -c 2 s V5p + ( Vv + 4r ? /3)V 2 u L + f L , (6a) 

d t Sp + p V-u L =g, (6b) 



where we have used the identity V(V • u) = V 2 u L . After separating the force-density field into 
transverse and longitudinal components, f = f T + f L , the two flows become independent. The 
sound is associated with the longitudinal flow, which is coupled dynamically to the density — i.e., 
the sound modes characterize the propagation and relaxation of a density perturbation. 

The time evolution of a density perturbation is obtained by taking the divergence of Eq. (|6ap 
and then transforming it using the mass-conservation equation, Eq. (|6bj) . 

d 2 Sp = c 2 s V 2 5p + 2rV 2 3^p - b • V5(r)5(t). (7) 

In Eq. ([7]) r = 1 nv+v^v ^ anc j we h ave se t g = and the force density to an impulsive point force, 
f = b<5(r)(5(i). Note that the application of the impulsive point force generates a momentary 
mass dipole in the equation for the density, — (b<5(i)) • V<5(r). Equation ([7]) represents a damped 
wave equation with a damping parameter V. If T = 0, we retrieve a wave equation for 5p, 



and if c s = 0, we obtain a diffusion equation for dt5p. Upon Fourier-transforming the spatial 
coordinates, Eq. ([7]) becomes 



SfSpk + (2Tk?)dt6p k + (c 2 k 2 )5p k = -ik ■ b5(t), 



whose solution reads 



Sp k (t) 



-ik • he~ TkH 

^cjk 2 - r 2 k* 



sin (ty/c 2 s k 2 -T 2 k* 



(8) 



Ifc 2 

-2 



- r 2 fc 2 > 0, the density perturbation propagates as a decaying wave. If, on the other hand, 
c" s - Y 2 k 2 < 0, the sine becomes a hyperbolic sine, and the density perturbation propagates 
diffusively. The crossover between these two behaviors occurs above a wavelength ~ T/c s . 
For water T/c s ~ 1 nm, while a typical colloid diameter is a ~ 10 3 nm. Therefore, particle 
correlations arising from sound are mediated by essentially underdamped wave-like modes. In 
this large-wavelength limit, k~ l S> r/c s , Eq. ([5} further simplifies to 



-ik • be" rfc2 * 
5p k (t) « sin(c s kt). 



0) 



Figure [TJ shows a snapshot of the density perturbation, generated a certain given time after 
the application of an impulsive force inside an unbounded fluid. We compare the result of a 
lattice-Boltzmann simulation, where the force is applied to a rigid spherical particle, with Eq. (j^J) 
for a point force. The analytical calculation reproduces all the features of the simulation, with 
a quantitative discrepancy that can be ascribed to the finite particle size in the simulation. 
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Figure 1. Snapshot of the density perturbation generated by an impulsive force in an 
unbounded fluid. The force is applied at the origin in the positive x direction. Lattice-Boltzmann 
simulation result (where the force is applied to a rigid sphere of diameter a) is compared with 
the analytical one for a point force [Fourier-inverted Eq. (jS|)]. The snapshot is taken 0.074rr 
after the impulse, where rr = h 2 /V and h = 1.5a. 



2.2. Channel with slip boundary conditions 

For a fluid confined in a channel with slip boundary conditions, we neglect variations of the 
density over the channel cross-section and consider only the direction along the channel, x. 



Equation ([7J) is then rewritten as 



d 2 t 5p = c 2 s 5p" + 2T5d tP " - ^S'(x)5(t), 

where a prime indicates differentiation with respect to x, and an impulse is evenly distributed 
over the cross-sectional area of a square channel of side h. In the long wavelength limit the 
solution takes the form, 



Sp 



2c,h 2 



-(x-c a t) 2 /(irt) 



-(x+c s t) 2 /(4Tt)' 

v / 4vrrt 



(10) 



There are two signals propagating with the speed of sound in opposite directions. In addition, 
each signal spreads out diffusively with diffusion coefficient T. Figure [2] shows density 
perturbation snapshots generated by an impulse applied at the center and along a square channel 
with slip boundary conditions. Lattice-Boltzmann simulation results, where the impulse is 
applied to a spherical particle, agree well with Eq. (|10p . 
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Figure 2. Three snapshots of a density perturbation for a fluid in a square channel with slip 
boundary conditions. An impulsive force is applied along the channel axis in the positive x 
direction. Lattice-Boltzmann simulation results (where the force is applied to a rigid sphere of 
diameter er) are compared with the analytical ones for a force distributed uniformly over the 
square cross-section of side h = 1.5a [Eq. (jlOp ], The time unit is tt = h 2 /F. 



2.3. Channel with no-slip boundary conditions 

To account for no-slip boundary conditions at the channel walls, we introduce an effective friction 
term, (— po£u L ), to the right-hand side of Eq. (f6a|) . This leads to an extra term, (—^dfSp), in 
Eq. ([7]). Thus, neglecting again density variations transverse to the channel axis, we have the 
effective ID equation, 

d?8p = c 2 5p" + 2Yd t 5p" - £d t Sp - ^S'(x)6(t). 

h A 

For the friction coefficient (rate of momentum loss to the walls) we take £ = av/h 2 , where 
v = rj/po is the kinematic shear viscosity, and a a geometrical prefactor, dependent on the 



shape of the channel cross-section. (The value of a can be determined from the resistance of the 
channel to a steady pressure-driven flow, e.g., a ~ 28.454 for a square cross-section and a = 36 
for a circular one [HE].) 

In Fourier space the solution reads 



(—ik)e 



-(Tk 2 +av/2h 2 )t 



Sp k (t) = — / ^ S m(tJc 2 s k 2 -(Tk 2 + av/2h 2 y). (11) 

h h 2 ^c 2 k 2 -(Tk 2 + au/2h 2 ) 2 V V s V ' ' > K ' 

The corresponding signal propagates with underdamped oscillations only for a narrow range of 
small wavelengths, satisfying 

/ <rr~\ I I or \ 

(12) 

where D s = c 2 h 2 /(au) is the sound diffusion coefficient. Hence, the friction at the walls 
transforms the large-scale dynamics from oscillating to diffusive. In the large-wavelength limit 
Eq. ([TTj) reduces to 

/ b 

i / 




5 Pk (t)^(—)(-ik)e- D ° k2t , (13) 



which in real space reads, 



b 2 X e-^l^st 

Thus, at large distances the density perturbation evolves as a diffusive dipolar signal. 

Figure [3] shows snapshots of the density perturbation propagating in a square channel with 
no-slip boundary conditions. The wave-like propagation quickly disappears and the perturbation 
continues to evolve diffusively. Once again, the simplified ID theory is found to qualitatively 
reproduce all the features seen in the lattice-Boltzmann simulation |21j . 




Figure 3. Two snapshots of a density perturbation for a fluid in a square channel with no-slip 
boundary conditions. An impulsive force is applied along the channel axis in the positive x 
direction. Lattice-Boltzmann simulation results (where the force is applied to a rigid sphere of 
diameter a) are compared with the analytical ones for a force distributed uniformly over the 
square cross-section of side h = 1.5a (Fourier-inverted Eq. ([11]) . multiplied by a factor of 2 [2T]). 
The time unit is rp = h 2 /T. 



3. Velocity cross-correlation function 

Section [2] has been concerned with the dynamics of density perturbations in a viscous fluid. We 
now proceed to examine the effect of this dynamic response on particles suspended in the fluid. 

Hydrodynamics represents the fluid as a continuous medium. Its discrete structure is revealed 
when a large particle immersed in the fluid undergoes Brownian motion due to large number of 
collisions with the much smaller fluid particles. Fluid particles impart momentum to the large 
particle, and the particle, in its turn, returns momentum to the fluid by setting off large-scale 
fluid flows. These flows lead to fluid-mediated correlations between the large particles. The flows 
can be decomposed, according to Eqs. ([S]) and into transverse and longitudinal components. 
The present work is concerned with correlations caused by density perturbations, and so we 
focus on the role of longitudinal flows. As noted in Sec. [TJ transverse flows in a rigid channel 
are screened at distances larger than the channel width and do not contribute to long-range 
correlations. 

To visualize these flow-mediated correlations, we show in Fig. [5] velocity cross-correlation 
functions between two isolated particles, as obtained from lattice-Boltzmann simulations for 
three different systems: an unbounded fluid, a fluid in a square channel with slip boundary 
conditions, and a fluid in a square channel with no-slip boundary conditions. The cross- 
correlation function, 

c(t,d) = (y 1 (o)v 2 (t)) d , 

measures ensemble-averaged correlations between the velocities of two particles, V\ and V 2 , 
separated by a distance d. In all three cases we take the two velocities to be oriented along 
the line that connects the centers of the two particles. To avoid the time-consuming fluctuating 
simulations and averaging procedure, we employ instead a deterministic measurement of C(t,d) 
based on the linear response of the two-particle system. The equivalence of the two measurements 
is guaranteed by the fluctuation-dissipation theorem [22] . The procedure goes as follows. For 
t < both particles are at rest, Vi(t) = V 2 {t) = 0. At t = we apply an instantaneous force to 
particle 1 by assigning it a finite initial velocity, V\{t = 0) = Vo- As time evolves, V 2 {t) responds 
to the perturbation set off by particle 1. The velocity cross-correlation function is then obtained 
as 

C(t,d) _ 1 V 2 (t) 
k B T MVi($y 

where M is the particle mass and k^T the thermal energy. 

In Fig. 2] we see that in all three systems there is a finite incipient time required for the signal 
to reach particle 2. This time is roughly d/c s . The leading signal in the channels is larger than 
that in the unbounded fluid, because the perturbation is guided along the channel rather than 
spreading in 3D. At larger inter-particle distances, however, that signal is suppressed in the case 
of a rigid channel with no-slip boundary conditions because of the diffusive nature of sound. 
(See Fig.EJ) 

The correlations at long times are very different for the three systems. In the unbounded 
case the correlation is governed by a positive algebraic tail, ~ t -3 / 2 . This is a manifestation of 
the well known long-time tails arising from the 3D diffusion of vorticity (monopolar transverse 
flow) |22] , In the rigid (no-slip) channel correlations due to transverse flows are screened, and 
the long-time correlation is governed by a negative algebraic tail, ~ (— i -3 / 2 ). This is a result 
of the ID diffusion of sound (dipolar longitudinal flow) presented in Sec. [2j By contrast, in the 
case of a channel with slip boundary conditions, the longitudinal signal propagates acoustically, 
making the correlation vanish exponentially fast once the signal has passed particle 2. 

We next focus on the sound-mediated correlations in a rigid channel with no-slip boundary 
conditions. Figure [5] shows C(t,d) for different separations d. The structure is similar for the 
different separations. There is the incipient time required for the leading signal to reach the other 
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Figure 4. Velocity cross-correlation functions for two particles, as obtained from lattice- 
Boltzmann simulations of three systems: an unbounded fluid, a fluid in a square channel (side h) 
with slip boundary conditions, and a fluid in a square channel with no-slip boundary conditions. 
The particles have diameter a and mass M and they are separated by a distance d = 2a = 1.33/t. 
The time unit is rp = h 2 /T. 

particle. This signal is suppressed with increased separation because of sound diffusion. It is 
followed by a long-time negative algebraic tail, whose magnitude does not show any dependence 
on separation. On the one hand, this indicates a correlation of unrestricted spatial range. On the 
other hand, the backflow associated with the negative tail causes the particle to return almost 
exactly to its initial position; the time-integrated correlation (which is equal to the coupling 
diffusion coefficient of the particle pair) decays exponentially with d/h. This is how the two 
phenomena — long-range temporal correlations and short-range steady-state correlations — are 
reconciled. 

To quantitatively account for the long-time algebraic behavior of the correlations we need to 
relate the diffusive density perturbation, given by Eq. (|14p . to the flow velocity. This can be 
done either via the continuity equation, dt5p = —pod x u, or by applying Fick's law to the sound 
diffusion, pqu = —D s d x 5p. Both methods yield the same long-time flow velocity, 



Equating it with the velocity advecting particle 2 and applying the fluctuation-dissipation 
theorem, as explained above, we obtain the long-time velocity cross-correlation function [23], 



which has no dependence on the inter-particle separation d. The negative sign of the long-time 
tail is traced back to the dipolar shape of the diffusive density perturbation, Eq. (|T4"]h 

4. Collective dynamics 

Finally, we study the effect of sound on the collective dynamics of many particles. The discussion 
in this section is restricted to the qlD geometry of a rigid channel with no-slip boundary 
conditions. 
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Figure 5. Velocity cross-correlation functions for two particles in a rigid channel (no-slip 
boundary conditions), as obtained from lattice-Boltzmann simulations. The particles, situated 
on the channel axis, have diameter a and mass M. The channel has a square cross-section of 
side h = 1.5o". Correlations are shown for three different inter-particle separations. The time 
unit is rr = h 2 /T. The inset shows the absolute value of the correlations on a logarithmic scale, 
demonstrating the (—3/2) power law at long times, with a distance-independent amplitude. 

We consider the wavenumber-dependent particle current correlation function [22], 



where ./V is the number of particles and i,j are particle labels. Self-contributions to the sum, 
with i = j, are related to the velocity autocorrelation function of single particles, (V(0)V(t)}. 
As our interest here lies solely in cross-correlations, we subtract this contribution and consider 



To obtain J c from lattice-Boltzmann simulations we assume time-scale separation between 
the relaxation of the fluid and that of particle configurations. We generate configurations 
X N of N particles using a Monte Carlo scheme. For each configuration we assign to each 
particle an initial velocity from a Maxwell-Boltzmann distribution and then follow the time 
evolution of every particle in a lattice-Boltzmann simulation. This allows us to calculate, 
for each configuration, the velocity cross-correlation functions between all pairs of particles, 
(Vi(0)Vj(t)){X N }, as described in Sec. [3l The assumption of time-scale separation simplifies the 
expression for J c to 



where g(x) is the equilibrium pair correlation function of the particles, and <fi = Na/L (L 
being the channel length) is their linear fraction. The angular brackets (. . . )fast 

denote a 

double averaging procedure: first, for a given configuration, an average is taken over time 
that is sufficiently long for the fluid to relax but sufficiently short for the configuration to be 
considered "frozen" ; then another averaging is performed over many particle configurations with 
the constraint x\ — X2 = x. 

In Fig. [6] we plot the time evolution of J c (k,t), thus obtained from the lattice-Boltzmann 
simulations, for various wavenumbers k. As t — > 0, J c vanishes, as the inter-particle correlations 




(16) 



J C (M) = J(k,t) - (V(0)V(t)). 




in the compressible fluid require a finite time to develop. The presence of oscillations for 
small wavelengths, comparable to the particle size and channel width, reflects the wave- 
like propagation of sound over these small length scales. The oscillations vanish at longer 
wavelenghs, indicating the onset of diffusive sound propagation. At long times, the particle 
current correlation J c decays algebraically as i~ 3 / 2 . 
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Figure 6. Particle current correlation functions obtained from lattice-Boltzmann simulations 
for a suspension of particles (diameter a, linear fraction <j> = 0.72) in a square channel (side 
h = 1.5<r). The different curves, corresponding to different wavenumber k, are shifted vertically 
by 0.15 from each other, for clarity. The time unit is rr = h 2 /T. 



To further clarify the results of Fig. [6] we resort to a simple model, where g{x) is replaced by 
its dilute limit: g(x) = for \x\ < a, and g(x) = 1 for |x| > a. For the velocity cross-correlation 
function (Vi(0)V2(f))f as t{a;} we use the ID description of Sec.O From Eq. (fTTj) and the continuity 
equation, dfSpk = —ikpou^, we obtain the flow velocity u^. The fluctuation-dissipation theorem 
is then used to relate the flow to the velocity cross-correlations. This model yields 



J c (k,t) 
k B T ~ ba 



Uk{t) - I dxe lkx u(x,t) 



(17) 



The qualitative change in J c as its oscillations vanish at long wavelengths is reproduced by the 
first term. The onset of the overdamped diffusive regime for small wavenumbers, lying outside 
the range delineated in Eq. (|12p . should occur for the parameters used in the simulation at 
ka < 0.18, which is in line with Fig. [6l The second term in Eq. (|17[) supplies the function J c 
with the algebraic decay ~ t _3//2 , which at long wavelength is positive. 

Another function that characterizes the collective dynamics of particles is the dynamic 
structure factor |22l [M] , 

1 N 

S (k,t) = -Lj2( eik(xm ~ Xj{t)) )- ( 18 ) 

At times much shorter than the configurational relaxation time it decays exponentially with 
time and can be written as 

S(k,t) = S (k)e- tk2H ^ s ^ k \ 



where So(k) = S(k, t = 0) is the static structure factor. The hydro dynamic factor H(k) accounts 
for the effect of fluid-mediated correlations on the dynamic structure and is given by [23] 



1 N I 

H(k) = Aelf + J^Y,\ D iji XN } eik(Xl 



-Xj) 



where Dij{X N } are the configuration-dependent pair diffusion coefficients, and -D S eif = (Da) 
is the self-diffusion coefficient. Note that H(k) depends only on steady-state properties of the 
fluid and, hence, is unaffected by sound. 
Recalling the Green-Kubo relation [22], 



dt ( Vi(0)Vj(t) ){x} 



Dij(x) 



and assuming time-scale separation, we relate H{k) to the current correlation function, 



H(k) = / dtJ(k,t). 
Jo 



To examine the effect of sound we therefore use the functions 
H(k,t) = 



t r-t 

dt'J(k,t'), H c (k,t)= / dt'J c (k,t'), 
o Jo 

such that in the limit t — > oo H(k, t) reduces to H(k), and H c (k, t) to H(k) — D S eif [22]. Figure[7| 
shows the dependence of H (k, t) on k at various times. As t increases, H(k, t) converges to its 
limiting form, H(k), and eventually all sound effects vanish. However, at any finite time there 
appears a feature that is absent in unbounded suspensions — a sharp peak at k — > 0. 
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Figure 7. Temporal hydrodynamic factor as a function of wavenumber at various times. The 
system is the same as in Fig. [6j 



To see that the small k feature is due to long-range sound-mediated correlations, we return 
to the simplified ID description. Using Eq. (fT7|) and the continuity equation, dtSpk = —ikpQiik, 
we obtain 



H c (k,t) 
k B T 



bapo 



(-ik) 



1 r dkl sm{cj{k' -k))5p kl {t) 



7T 



k'-k 



-ik') 



(19) 



where 8pt is given by Eq. (jlip [or by Eq. (|13j) for small k\. At t — > oo this expression vanishes, 
indicating that within the simplified ID description all steady correlations between particles 
confined in the channel disappear, and H(k) = D SD \{. However, at any finite t, we have a range 
of small wavenumbers, k 2 < (D s t) , for which Eq. ()19p reduces to 



Dq a \ \firD s t 

where Dq = ksT '/ \3Trrja) is the Stokes-Einstein self-diffusion coefficient in an unbounded 
fluid. Equation (|20p contains two interesting terms originating from sound diffusion. The first 
contributes at k — > a time-independent constant. The second adds a k-independent negative 
long-time algebraic tail, ~ (— 1~ 1 ^ 2 ), which derives from the time-integrated i -3 / 2 tail of the 
current correlation function J c (k,t). 

In Fig. [8] we replot the function H(k,t), focusing on the small k region, and compare it with 

the fit = ^i-e~ Dak2t |21j . The agreement between theory and simulation is satisfactory, 

both exhibiting the small-A; peak. 




Figure 8. Temporal hydrodynamic factor for small wavenumber at various times. The system 
is the same as in Fig. [6j Data points represent lattice-Boltzmann simulation results. Solid lines 
show the expression — = ®ink e -D s k 2 t _ 



5. Conclusion 

In the present work we have tried to establish the role of fluid density perturbations (sound) 
in mediating dynamic correlations among particles in a channel at the pair and the collective 
levels. We have been motivated in this investigation by two considerations: (i) the hydrodynamic 
effects that are not associated with sound (due to momentum and steady pressure) are screened 
in a rigid channel and have no bearing on long-range particle correlations; (ii) the density 
perturbations in a channel do not decay within a well defined relaxation time (as in an unbounded 
fluid) but rather relax algebraically, thus giving sound a long memory. 

The effects of sound on particle correlations in a rigid channel have been found to be quite 
striking. At the two-particle level, the velocity cross-correlation function is independent of 
inter-particle distance (i.e., unrestricted in its range) and has a long-time negative algebraic tail 
~ (— t -3 / 2 ) (similar to the single-particle autocorrelation function [5l El [7j)- At the collective 
level of a large particle assembly (qlD suspension), the particle current correlation function has 



a similar long-time algebraic tail. We have characterized the collective correlations by a time- 
dependent hydrodynamic factor, H c (k,t), which converges to the ordinary hydrodynamic factor 
at infinite time [25] . The sound-mediated correlations in the channel make this function exhibit 
an anomalous peak at small wavenumbers and a long-time algebraic temporal decay ~ (— t" 1 / 2 ). 
Thus, in summary, the common picture of negligible, exponentially screened hydrodynamic 
interactions among particles in a rigid channel is somewhat misleading; strictly speaking, it 
applies only at steady state. 

At the same time, the predicted sound-mediated correlations are weak and may be measurable 
only over sufficiently short times. The correlation is inversely proportional to the speed of sound 
[see, e.g., Eq. (|15|) ]. which makes its amplitude small, and the sound diffusion is fast — for 
water in a micron-wide channel, D s = c^h 2 /{au) ~ 0.1 m 2 /s. For example, let us examine the 
correlation between the observable Brownian displacements of two distant particles along the 
channel, as characterized by their pair diffusivity, 

(AxiAx 2 ) _ 1 [* f A , ru , l j. _ k B T 



2t t J J c s hpo^7rais 



For a micron-wide channel filled with water, according to Eq. (|2ip . measuring D12 of order 
10~ 3 /iin 2 /s would require a temporal resolution of order 10 -7 s. Thus, the predicted sound- 
mediated effects might be observable using experimental techniques of high temporal and spatial 
resolution, such as diffusive wave spectroscopy [12] or fast optical tracking [14L I26j . The fact that 
-D12 is independent of inter-particle distance may greatly facilitate the acquisition of statistics 
in such fast-tracking experiments. 

The short-time relevance of sound-mediated correlations was noted in earlier works on 
unconfined suspensions \12\ IT3j I14j . We end by underlining the key differences between those 
correlations and the ones in a rigid channel addressed here. First, whereas in an unconfined 
suspension sound propagates as an underdamped wave until it is scattered by a sufficient 
number of particles [T2], in a channel, because of scattering from the boundaries, sound becomes 
diffusive as soon as the distance gets much larger than the channel width. Second, in an 
unconfined suspension sound- mediated correlations are a sub-dominant effect, quickly taken over 
by vorticity diffusion (transverse modes) [12^ I13j. By contrast, in a rigid channel correlations 
due to transverse flow are cut off, leaving longitudinal flow as the sole mechanism of long-range 
correlations. Thus, any correlation which could be resolved at inter-particle distances much 
larger than the channel width should be associated with sound-mediated effects. 



Acknowledgments 

This research has been supported by the Israel Science Foundation (Grants No. 588/06 and No. 
8/10). 



References and Notes 

[1] Happel J and Brenner H 1983 Low Reynolds Number Hydrodynamics (The Hague: Martinus Nijhoff) 

[2] Cui B, Diamant H and Lin B 2002 Phys. Rev. Lett. 89 188302 

[3] Brinkman H C 1949 Appl. Sci. Res. A 1 27 

[4] Ramaswamy S and Mazenko G F 1982 Phys. Rev. A 26 1735 

[5] Hagen M H J, Pagonabarraga I, Lowe C P and Frenkel D 1997 Phys. Rev. Lett. 78 3785 
[6] Pagonabarraga I, Hagen M H J, Lowe C P and Frenkel D 1998 Phys. Rev. E 58 7228 
[7] Pagonabarraga I, Hagen M H J, Lowe C P and Frenkel D 1999 Phys. Rev. E 59 4458 
[8] Frydel D and Rice S A 2007 Phys. Rev. E 76 061404 
[9] Felderhof B U 2011 J. Chem. Phys. 134 024505 

[10] Cui B, Diamant H, Lin B and Rice S A 2004 Phys. Rev. Lett. 92 258301 

[11] Diamant H 2009 J. Phys. Soc. Jpn. 78 041002 

[12] Ladd A J C, Gang H, Zhu J X and Weitz D A 1995 Phys. Rev. E 52 6550 



[13] Bakker A F and Lowe C P 2002 J. Chem. Phys. 116 5867 

[14] Henderson S, Mitchell S and Bartlett P 2002 Phys. Rev. Lett. 88 088302 

[15] Felderhof B U 2010 J. Fluid Mech. 644 97 

[16] Frydel D and Diamant H 2010 Phys. Rev. Lett. 104 248302 

[17] Navardi S and Bhattacharya S 2010 J. Math. Phys. 51 043102 

[18] Felderhof B U 2010 Physica A 389 5602 

[19] Felderhof B U 2011 J. Fluid Mech. 668 100 

[20] Landau L D and Lifshitz E M 1987 Fluid Mechanics 2nd edition (Pergamon: Oxford) 

[21] The fields (density, flow velocity) in the ID model correspond to cross-sectional averages of the actual fields, 
whereas in the simulation we measure quantities at the center of the channel. This makes the simulation 
result larger than the analytical one by a factor of « 2. (For a parabolic velocity profile in a circular tube 
the cross-sectional average is exactly 1/2 the central value [9].) 

[22] Hansen J P and McDonald I R 1986 Theory of Simple Liquids 2nd edition (London: Elsevier) 

[23] Because of cross-sectional averaging [5T], this result of the ID model is smaller than the corresponding long- 
time behavior in the simulations by a factor of ~ 4. The additional factor of ~ 2 arises from the continuity 
equation (or Fick's law), relating the cross-sectional averages of density (a roughly uniform transverse 
profile) and velocity (a roughly parabolic profile). 

[24] Pusey P N 1991 in Liquids, Freezing and Glass Transition IL eds J P Hansen, D Levesque and J Zinn- Justin 
(Amsterdam: Elsevier) pp 763-942 

[25] A time-dependent hydrodynamic factor was introduced before, using a different definition [T2]. The relation 
between our function, Hc 2 \ and the one used in Ref. [12], is H^\k,t) = t" 1 f* dt' (k,t'). 

[26] Huang R, Chavez I, Taute K M, Lukic B, Jeney S, Raizen M G and Florin E-L 2011 Nat. Phys. 7 576 



